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In this work we investigate the chiral symmetry breaking in the o-deformed 

00 ; 

I version of the NJL model and its consequent mass generation mechanism. We show 

. that the deformation of the NJL model, in the mean field approximation, may take 
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into account correlations that go beyond the mean field and, in a certain limit, 



^ ■ approaches the more realistic lattice calculations. 
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Quantum algebras, or simply q-deformed algebras, known in the case where just a 
single deformation parameter is introduced, have received some attention last years in the 
study of many-body problems. Because they can provide us with a class of symmetries 
that is richer than the usual Lie algebras, they are appropriate to describe physical systems 
and (or) models which are not properly described by the last ones. On the other hand, 
the introduction of the g-parameter in the theory can be viewed as a way to take into 
account correlations in many-body systems. As examples, we mention a set of works p], 
related with the deformation of the su(2) algebra (also called su g (2) algebra) in a time 
dependent Hartree-Fock (TDHF) approximation to solve the Lipkin model. Following 
similar reasonings, the quon algebra, which describes particles whose statistics interpolates 
between the bosonic and fermionic ones 0, was then used to perform a g-deformed boson 
analysis of the random-phase-approximation (RPA) solution for the Lipkin and the two- 
level pairing models M. In that case there is a strong indication that the introduction 
of correlations through the deformation of the RPA bosons can also be interpreted as 
a restoration of a broken symmetry caused by using the regular (approximated) RPA 
solution. This is manifested in the smearing out of the phase transition at the so-called 
critical point pi. 
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Recently, the dynamical breaking of chiral symmetry in the Nambu-Jona-Lasinio 
(NJL) model || and the corresponding generation of a dynamical mass for quarks, were 
analyzed at the light of a g-deformed version of that model f7j. Once g-deformation ap- 
pears as a powerful tool to study symmetries and their generalizations, this seems to be 
a quite appealing investigation. For that purpose, a particular prescription for the defor- 
mation was chosen, based on reference || and two important conclusions of that analysis 
were drawn: the phase transition of the NJL model continues to be sharp and occur at 
the same critical interaction point, when compared with the non-deformed case, which 
also means that any possible explicit chiral symmetry breaking terms in the Lagrangian 
cannot be simulated by the deformation procedure used; the main modification caused by 
the deformation is an enhancement of the condensate, i.e., the quark mass is effectively 
increased for a given value of the interaction constant. 

In the present work, it is our intention to proceed further with the analysis introduced 
in reference [7). We propose a method to solve the NJL model based on a Hamiltonian 
density written in terms of the su(2) operators. As we show bellow, this is in fact possible 
for the NJL model and corresponds to a J=l angular momentum algebra ||, once we are 
restricted to the well-known BCS type ansatz for the variational solution in the Hartree 
approximation. With that result in hand, it is then straightforward to build an su g (2) 
version of the NJL model. This has the advantage that it is then possible to exploit 
different g-deformation schemes [TO]. Moreover, this procedure allows us to generalize 
some previously obtained results and to explore all features introduced by g-deforming 
the model. 

We start with the Lagrangian density of the NJL model |j) 

L = iTid^ + G ((^) 2 + (feWO 2 ) , (1) 

where ip represents the quark field amplitude, the second term represents the interaction 
mechanism and G is the coupling constant. From (TJ), the corresponding Hamiltonian 
density can be obtained and reads: 

U = -ifj.Vip - G ((^V) 2 - (?i7s^) 2 ) • (2) 

As it is well known, a solution for the above problem can be found in the Hartree 
approximation, using a BSC-like variational state of the form: 

\NJL >= H (cosfl(p) - ssin^(p)6 t (p,s)rf t (-p,s)) |0 >, (3) 
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where p is the vector momentum, s =±1 are the helicity state eigenvalues and the particle 
and anti-particle operators are defined such that 

6(p,s)|0>=0 , d(p,s)|0>=0. 

In terms of the creation and annihilation operators, we expand the field operators at 
t=0, yielding 

^(x, 0) = E / (S3 (6(P, S MP, s)e ip ~ + ^(p, s)v(p, s)e-^) , 

where u(p, s) and v(p, s) are the normalized spinor eigenfunctions for particles and anti- 
particles with momentum p and helicity s. 

At this point, we introduce the following angular momentum operators: 

J +jP = V24 , J_ tP = V2A p , J 0p = ^P-l, (4) 

where 

4 = 7f ( &t (P' +) rft ("P> +) - &t (P, ~)d\-p, -)], (5) 

= E[ fot (P. S )KP, s) + dt(p, s)d(p, s)\, (6) 

s 

the operators </+, p , </-, p and Jo,p obey the usual su(2) commutation relations and 

(J+, P ) n |0 >= 0, if n > 2, J_ p |0 >= 0, J ,p|0 >= -|0 > . (7) 
The variational ansatz defined in equation (|3]) can then be written as: 



\NJL >= Afe * |0>, (8) 



with 



A^= FTT5-, (9) 

n(i + ¥) 

p 

where £ p = a/2 tan ^(p). In terms of those operators, the hamiltonian (Q) (in the Hartree 
approximation) can be put in the form: 

H = J d 3 xH = J -0^ {p 2 J , p + M (J +)P + J_, p )} , (10) 

where M represents the mass of the arising condensate. 



We now minimize the mean value of the Hamiltonian with the above NJL state, using 
6(p) as the variational function and end up with the well known NJL solution for the 
condensate, i.e., < ijnj) >= —M/(2G). At this point we introduce the deformation in the 
model. Once the algebraic su(2) structure above underlies the model, there exists a direct 
procedure to extend the solution to the su g (2) formalism ]1(J. The g-deformed angular 



momentum operators obey the following commutation relations: 

[ Jo, p , J ±lP ] = ±J ±)P , [J +iP ,J-, p ] = (gr) J - Jo -[2Jo, p ], (11) 



where J = J J\ v + Jl p + J$ p and [X] represents the deformed version of the operator 
(or c-number) X. Its definition is not unique. Throughout this paper we use the following 
definition pi: 

X —Y 

m = q - — — , (12) 

q — r 1 

and only cases where just one deformation parameter is introduced is considered, namely, 
r = q and r = 1, respectively. Note that in the first case (r = q), q can be either a real or a 
complex number (q = e lT , r being real), and in the second one (r = 1) it must be real. To 
proceed with the deformed solution we introduce a g-deformed variational ansatz, which 
is: 



\NJL > g =Af q e q p |0 >,) (13) 



where 



Mq n(i + ^ + (¥) 2 ) (14) 

and the g-exponential is defined as e q {ax) = or^™ > with [n]! = [n][n — 1] [2][1] . 

n=0 

Again, we minimize the mean value of the Hamiltonian flIU|), using the variational 



state defined in (|T3|) and bearing in mind that the angular momentum operators obey the 
deformed commutation rules. After a straightforward calculation we obtain the modified 
gap equation: 

Af _ irr oi M 3 P f tan%)(l + tan 2 %)) \ 

[ l J (27r) 3 \[l + [2]tan 2 ^(p)+tan 4 ^(p)]J ' { } 

where 6{p) is a solution of the equation 

tan 6 d{p) + ^ tan 5 9{p) + (3 - [2]) tan 4 9{p) + ^ tan 3 6(p) 
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-(3- [2]) tan 2 6(p) + ^ tan 9{p) -1 =0. 

At this point it is useful to make some comments about the above equations. Our 
choice for a g-deformed exponential in the NJL ansatz is somewhat arbitrary, in the sense 
that we could have introduced the deformation just by modifying the su(2) commutation 
rules. If we had chosen to write the variational state (|T3"D in terms of a regular exponential 
function, our results would reduce to: 

M * _ iG * f d 3 p I M* \ 

J (2VT)3 I 0,2 + (M *)2 J ' 1 ) 

where M* = \f^M and G* = G[2). As [2] -> 2 when q -> 1, this last result also 
makes clear that we reproduce the right limit for the gap equation [11]. We have verified 
numerically that the difference between the results obtained from the gap equations (15) 
and (16) is negligeable, unless we investigate the behavior of the dynamical mass for q 
values very far from 1, which is not the case in the present study. Thus, the introduction 
of the deformed exponential, instead of the regular one, in the variational ansatz, does 
not affect the results discussed next. 

In figures la and lb we show the condensate as a function of the inverse coupling 
constant for the two kinds of deformation considered here, i.e., r = q and r = 1 respec- 
tively. In the first case, q can acquire any complex values, but in the second one it must 
be real. We have chosen two values in each case in order to show the main modifications 
introduced by the deformation. In figure la, for r = q, we see that for q real we obtain 
an increase of the condensate, in agreement to what was found in |7|], and for a complex 
q the condensate decreases, in relation to the condensate obtained for q — 1. This is not 



surprising, once according to |L5| and [|16|] , complex q- values mimic a repulsive interaction 
while real ones play the role of an attractive interaction. A similar effect can be seen when 
r = 1, for which only real g-values are allowed. The results are displayed in figure lb. In 
that case, however, for q < 1 a repulsive effect is simulated, while for q > 1 an additional 
attraction appears. In all cases we notice that the critical value for the coupling constant 
depends on the deformation parameter. It can be shown that the phase-transition point 
is given by G q c = with G c = where A is the well-known NJL cutoff. Throughout 
this paper we have used A = 600MeV. This is a feature not present in the result obtained 
in j7| and we return to this point latter. 

It is also worth mentioning that there is a certain degree of ambiguity in the way a 
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physical system can be g-deformed. This problem has already been extensively discussed 



in the literature |T3|, [14]]. In this way, equation fllPf) can be substituted by 

H = J d 3 xH = J -j^- {p [2 J 0jP ] + M (J +)P + J_, p )} . (17) 

before the minimization procedure is performed within the deformation scheme. As the 
only modification refers to the way we write the kinetic term, i.e., we replace 2 Jo by [2 Jo], 
the corresponding results are called the deformed kinetic (DK) results and the previous 
results obtained by minimization of the Hamiltonian ( JlQ| ) are called non-deformed kinetic 
results (NDK). A comparison between the two approaches is shown in figures 2a and 2b, 
for r = q and r = 1 respectively. In figure 2a, we can see that the DK results simple turn 
off de q dependence in G c , yielding the same behavior obtained in reference 0. However, 
for r = 1 the DK results can lead to strong quantitative modifications, specially in what 
concerns the phase transition point, as is made clear for q = 2 in figure 2b. 

In order to get some insight on the physical background embedded in the above results, 
we have decided to make a qualitative comparison with more elaborate calculations in the 



NJL model |17]]. In that work, a lattice Monte Carlo simulation and a Schwinger- Dyson 
(S-D) calculation are applied to obtain the condensate and compared with the usual gap 
equation solution. The first important result that emerges is the fact that the value of 
the critical coupling constant ( G c ) depends on the approximation used. The better the 
approximation, the bigger the critical coupling constant. Moreover, both the Monte Carlo 
simulation and the S-D method gives a smaller value for the condensate for any given value 
of the coupling constant, compared with the gap equation result, as can be seen in figure 
3 of |I7[. These two features can be simultaneously reproduced by our g-deformed NJL 
calculation for a complex q- value (r = q) or for q < 1 (r = 1). In figures 3a and 3b we 



plot the condensate as a function of for a better comparison with the results in (17 



In summary, we have introduced deformation in the NJL model by means of the angu- 
lar momentum operators, which are the generators of the su q (2) algebra. The condensate 
was then obtained by minimizing the deformed Hamiltonian with the help of a BCS ansatz 
for the q variational state. The results show that the phase transition is never suppressed, 
but the behavior of the condensate depends on the definition of the deformed quantity and 
on the numerical value of q. Two different prescriptions to deform the Hamiltonian and 
two different ways of introducing the g-deformed variational state were used. A qualita- 
tive comparison was made with calculations that go beyond the familiar Hartree solution 
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T7|. We can obtain the same behavior as these sophisticated calculations for any q- value 
which simulates a decrease in the residual attractive interaction. Our calculations give 
rise to a q dependence in the critical strength G c . This dependence is consistent with 
the repulsive character that can be mimiced by the deformation and goes to the right 
direction as compared with realistic calculations for the condensate. 
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fig. la. The condensate as a function of G c /G,for r=q in the deformed quantity definition. 
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fig2a. Same as fig. la for the deformed kinetic (KD) operator 
and for the non deformed kinetic (NDK) operator. 
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fig.2b.Same as fig. 2a with r=l. 
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fig3a. The condensate as a function of G/G for r=q. 



